Scattering lengths and universality in superdiffusive Levy materials 
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We study the effects of scattering lengths on Levy walks in quenched one-dimensional random and 
fractal quasi-lattices, with scatterers spaced according to a long-tailed distribution. By analyzing 
the scaling properties of the random-walk probability distribution, we show that the effect of the 
varying scattering length can be reabsorbed in the multiplicative coefficient of the scaling length. 
This leads to a superscaling behavior, where the dynamical exponents and also the scaling functions 
do not depend on the value of the scattering length. Within the scaling framework, we obtain an 
exact expression for the multiplicative coefficient as a function of the scattering length both in the 
annealed and in the quenched random and fractal cases. Our analytic results are compared with 
numerical simulations, with excellent agreement, and are supposed to hold also in higher dimensions. 
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I. INTRODUCTION 

Diffusion in heterogeneous and porous materials [ljj, 
composed of two or more types of regions with very dif- 
ferent diffusion properties, can be described as a sequence 
of independent scattering events occurring in the hard- 
scattering part of the material, followed by long jumps 
performed at almost constant velocity in non-scattering 
regions. The single-particle dynamics thus amount to a 
random walk, where each step length I is a random vari- 
able with a given probability distribution A (7). As it is 
well known, if I has a finite variance, the ensuing pro- 
cess follows the standard laws of Brownian motion. If 
instead the material is very heterogeneous on all scales, 
the step length distribution may become heavy-tailed and 
when (I 2 ) diverges the diffusion can become anomalous 
0. H| ■ Whenever the detailed structure of the underlying 
scattering media can be ignored (annealed disorder), the 
dynamics is only ruled by the behavior of X(l) for large I 
and the model corresponds to a standard Levy walk [J- 
6]. If instead one takes into account that the steps are 
correlated by their mutual positions in the sample, the 
step length distribution represents a quenched disorder. 
It is argued that quenched and annealed disorder differ 
in many respects |7Hlfj|. 

A particularly interesting application is multiple scat- 
tering of light in disordered media, which can often be de- 
scribed as a random walk process, analogous to the Brow- 
nian motion of massive particles [111 ] . In that case, if scat- 
tering elements are homogeneously distributed in space, 
one usually observes diffusive transport. However, when 
the local density of scattering elements is strongly inho- 
mogeneous, the optical transport properties of a material 
can change dramatically and lead to superdiffusion, as ex- 
emplified for instance by photon diffusion through clouds 
1 121 . More recently, artificial materials have been assem- 
bled in the lab that allowed for an unprecedented ex- 
perimental demonstration of light superdiffusion (l3l . [l4| . 



Interestingly, in experiments a change in the optical den- 
sity of the diffusive media allows to tune the so called 
scattering length (or scattering mean free path) [TJ, [ll| , 
which is a measure of the probability of experiencing 

Also Monte Carlo simulations 
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a scattering event 

[l5| indicate that this is an important parameter, that 
may significantly affect the observability of asymptotic 
scaling regimes. It is thus important to asses the role 
of the scattering lengths on scaling properties of rele- 
vant observables, starting with probability distributions, 
to interpret correctly experiments and simulations in in- 
homogeneous scattering media. While in pure diffusive 
processes a varying scattering length can be simply en- 
coded in a trivial change of the variance of the Gaussian 
distributions, superdiffusive process have not been inves- 
tigated in details. 

Many relevant features of experiments on scattering in 
inhomogeneous media can be described as a random walk 
in a quenched, long-range correlated environment. The 
simplest case consists of a free particle moving through a 
one-dimensional array of barriers whose spacing is power- 
law distributed |16l42fj |. To address the question dis- 
cussed above, in this work we extend such models to 
the case in which the transmittance through each bar- 
rier can be tuned to be different from 1/2. This models a 
situation in which the velocity in the scattering media is 
not fully randomized at each collision. Changing the bar- 
rier transmittance is thus akin to changing the scattering 
length in the diffusive portion of the material. 

In the following, we will show that the probability dis- 
tribution for the random walker to be a time t in the 
point r exhibits a scaling form and that the effect of 
varying the transmittance can be reabsorbed in the mul- 
tiplicative coefficient of the scaling length of the process. 
Numerical simulations confirm this superscaling behav- 
ior and evidence that not only the dynamical exponents 
are universal, as expected from the general scaling frame- 
work, but also the scaling functions are unchanged. No- 
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tice that the latter are in general not Gaussian in these 
systems. Moreover, we obtain an analytic expression for 
the multiplicative coefficient of the scaling length as a 
function of the transmittance both in the annealed and 
in the quenched models; in the annealed case the an- 
alytic form is calculated by solving directly the master 
equation of the process; in the quenched models it can 
be derived, within the scaling hypothesis, by consider- 
ing an extended Einstein relation between the stationary 
conduction and the (super)diffusion. Interestingly such 
coefficient turns out to be different in the quenche and 
in the annealed models and, in the latter case, it is in- 
dependent of the transmittance when superdiffusion is 
present. This evidences another relevant difference be- 
tween the two approaches. Finally, we apply the same 
scaling picture, as a function of the transmittance, to the 
time resolved transmitted intensity through a finite sam- 
ple of length L and to the total transmission, which are 
the quantities usually measured in experiments. 

The paper is organized as follows: in the next section 
we define the random and fractal quasi-lattices models 
with arbitrary transmittance, introducing the scaling pic- 
ture and the superuniversal behavior. In Section IIHI we 
calculate the exact expression, as a function of the trans- 
mittance, of the multiplicative coefficient for the scaling 
length in the annealed case. In Section HVl we turn to the 
quenched random and fractal quasi-lattices and we ob- 
tain an exact expression for the coefficients in these cases. 
In Section [V] we analyze the superscaling framework for 
random structures with Levy parameter < a < I, a 
case with non Gaussian scaling functions. Then, we fo- 
cus on random structures with 1 < a < 2, where the 
Gaussian scaling function obeys superscaling, but a sub- 
leading term determines the moments of the distribution. 
We evidence that also this subleading term can be stud- 
ied in terms of a change in the multiplicative constant 
for the scaling length. We finally verify the superscal- 
ing framework in the deterministic case with < a < I 
where the scaling function presents log-periodic oscilla- 
tions [l9|,[2(|. In Section lVTl we apply our superscaling ap- 
proach to a different physical quantity, the time resolved 
transmission intensity, which is relevant for a realistic ex- 
perimental setup in inhomogeneous optical media. In the 
last Section we present our conclusions. 



II. LEVY WALKS ON RANDOM AND 
FRACTAL QUASI-LATTICES 

We consider a continuous time random walk on one- 
dimensional structures (see Fig. [TJ), where point-like scat- 
terers are spaced according to a Levy distribution. The 
walker moves with constant velocity v in between each 
two consecutive scatterers. When it arrives at a scat- 
terer, it can reverse its direction of motion (v — > —v) 
with probability (1 — e)/2 (with — 1 < e < I). Clearly, \v\ 
is conserved during the evolution and we can set \v\ = 1 
without loss of generality. 



We investigate two types of structures. The first is 
random and (upper panel of Fig. [1} the probability for 
two consecutive scatterers, labeled by the indexes j and 
j + 1, to be at distance r is [HI, [l8[ 

A(r) = ^-, re[r ,co), (I) 

where a > and ro is a cutoff fixing the scale length of 
the system. The second type is a class of deterministic 
quasi-lattices (lower panel of Fig. [T]), built by placing the 
scatterers on generalized Cantor sets [Til HO]- Each set, 
and the ensuing step length distribution, is defined by the 
two parameters n u and n r used in its recursive construc- 
tion. The former represents the growth of the longest 
step when the structures is increased by a generation, so 
that the longest step in a structure of generation G is 
proportional to n„ ; n r is the number of copies of gener- 
ation G — 1 that form the generation G, so that the total 
number of scatterers in the generation G is proportional 
to (see Ref. 19] for details). For this second type of 
structures, the role of the exponent a of the random case 
is played by a = logn r /logn„ 20}. 

In the random case, we will average over different re- 
alizations of the structure and we will consider averages 
taken over processes starting from scattering sites. For 
quenched Levy processes, it is known that different aver- 
aging procedures can lead to different behaviors. More- 
over, properties arising from averages taken over pro- 
cesses starting in any point are different @, [To], [H, [2(| ■ 
In the deterministic case, we consider averages performed 
over random trajectories starting from a given point, e.g 
the origin evidenced in Fig. [T] 

The main quantity we are interested in is the proba- 
bility for a walker to be at time t a distance r from the 
starting point, which we denote by P a . £ (r,t) to empha- 
size the dependence on the two basic parameters of this 
class of models, a and e. The process on a Levy structure 
with given a can be described by introducing a scaling 
function f a {x) and a scaling length £ £ (t) growing with 
time [llj]. The scaling form for P a e (r,t) can be written 
as: 

P a , E (r,i) =£- 1 (t)f a (r/i e (t)) + h a , s (r,t) (2) 

where h a ^ e (r,t) is a function that vanishes in probability 
for large times: 

r-vt 

lim / \P a>e (r,t)-t- 1 (t)f a {r/l E {t))\dr = (3) 

In [18J, the scaling form (J2j , <j3j) , has been analyzed in 
details and it has been shown to hold for e = 0. More- 
over, the growth of £o(t) has been shown to follow the 
asymptotic law £ (t) ~ for < a < 1 and 

£a{t) ~ t 1 / 2 for a > 1. The presence of the subleading 
function h a fi(r, i) and of long tails in f a (x), induce a non- 
trivial scaling of higher-order momenta, (r^jt)) £ p (t), 
leading to strongly anomalous diffusion [2l| . 



3 



A 



write 



2 3 



,(!+£) 



4 H 6 



B 



FIG. 1: (a) Example of a random structure, where scatterers 
are distributed according equation fTJ. A particle with veloc- 
ity v is transmitted with probability (1 + e)/2 and reflected 
with probability (1 — e)/2; (b) A deterministic Cantor like 
Levy quasi-lattice, characterized by n r — 2 and n u — 4. 



Here we will evidence that, within a superscaling 
framework, the dynamical exponents and the scaling 
functions are independent of the transmittance e. In par- 
ticular, fa(-) is independent of e and the scaling length 
reads: 




if < a < 1 
if 1 < a 



(4) 



Notice that i e (t) is defined up to an arbitrary multiplica- 
tive constant, that in numerical simulation will be fixed 
equal to one for e = 0. Within the hypothesis that f a {-) 
is independent of s, in the next sections we will obtain 
an analytic expression for A e both in the annealed and 
in the quenched case. 



III. TRANSMITTANCE AND ANNEALED 
LEVY WALKS 



Before discussing the effect of the transmittance e in 
quenched systems, we consider the annealed case, where 
the length of the ballistic stretches is chosen randomly 
from the distribution ([1]) independently at each scatter- 
ing event. In this case the topology of the system is not 
taken into account and only the distribution of the steps 
influences the dynamics. In this simpler situation, an 
analytical approach is feasible and the behavior of the 
scaling length as a function of the transmittance can be 
determined. Let us introduce P + (r,t) and P~(r,t) as 
the probabilities of being in r at time t arriving respec- 
tively from the left and from the right. In the persistent 
random walk approach, recalling that \v\ = 1, we can 



P+(r,t) 



1 + e 



P + (r -r't 



1 - £ 

H ^- p ~( r ~r ,t- r') 

p-(r,t) 
1-e 



l+e 



\(r')dr' + -6{r,0)S(t,0) 



P~(r + r',t-r')+ 



P + (r + r',t-r') 



X(r')dr' + -6(r,0)6(t,0). 



(5) 



By applying a Fourier transform in space and time in 
([5]), we obtain an expression for P(k, ui), that is the trans- 
formed of Pa,e{r, t) — P + {r 1 t) + P^ (r,t). In particular in 
the asymptotic regime for large space and times, i.e. for 
ui and k going to zero, we have the following estimates: 



if < a < 1 
if 1< a < 2 



f (Cic^ + Cafc")- 1 
P a , e (h,uj) = hc^ + C±k°<)- 1 

[(C 5 w + C 6 ±i^fc 2 ) 1 if2<a 

(«) 

where the complex constants C\ . . . CV depend on the dis- 
tribution A(r) (i.e. on a) but are independent of e. The 
asymptotic scaling form of P a . £ (r, t) can then be obtained 
by inverting the Fourier transform in ([5]). In particular 
we have: 



P a *(r,t)2il e {t)- 1 f a {r/l e (t)) 



with 




if < a < 1 
if 1 < a < 2 

if 2 < a 



(7) 



(8) 



where the function / Q (-) and proportionality constants in 
(|8|) are independent of e. 

As expected, the dynamical exponents in equation ([8jl 
coincides with the well known results for the annealed 
Levy walks with e — 0, established in [6]. In addi- 
tion, equations (|7I8[) show that both the dynamical expo- 
nents and the scaling function / Q (-) are independent of 
e, even in the non trivial case of anomalous diffusion i.e. 
a < 2. Moreover, quite surprisingly for a < 2, i.e. in the 
ballistic and superdiffusive cases, even the coefficient of 
the scaling length does not depend on the transmittance 
and therefore the whole asymptotic regime is indepen- 
dent of the value of e. On the other hand, for a > 2 
the same coefficient depends in a non trivial way both 
on e and on the step length distribution A(r). Indeed, 
C 7 = (2(r) 2 /(r 2 ) - 1) where (r fc ) = / r k X(r)dr. In par- 
ticular, for the step lengh distribution considered in ([T}, 
C 7 = {a 2 -2a- I) /{a - l) 2 . Notice that the diffusiv- 
ity diverges for e = 1 (perfect transmission), while for 
£ = — 1 it vanishes only for CV = 1, hence only when 
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FIG. 2: Value of the generalized diffusivity in ((8} as a func- 
tion of e. Symbols represent numerical simulations, while 
continuos lines are the analytical calculations. 



the step lengths do not fluctuate. Indeed, step length 
fluctuations induce diffusion even in the case of total re- 
flection. Fig. [^compares the analytical prediction for the 
coefficient in Equations ([5} with numerical simulations at 
different value of a showing an excellent agreement. 



master equation reads: 

PaM) = ^P-in+i) + ^P + (r fe+ i) (11) 

where Pa E ( r k) an d Pae( r k) represent the stationary 
probabilities of being at the scattering site r k arriving 
from the left and from the right, respectively. The solu- 
tion of equation (fTTj) is 



P+ e (r k )=ak-- +c 

1 — £ 



(12) 



where a and c are arbitrary constants. Imposing the par- 
ticle density at the borders P a , e (0) = P+ £ (0) + P+ e (0) = 
1 and P a . e {L) = P£ e (L) + P+ e ( L ) = O ^e get c'= 1/2 
and a — (2K)^ X where K is the number of scatterers 
between and L. Since the conductivity equals the total 
current flowing in the system we have 



C a<E (L) = P+ e (r k+1 ) - P~ e (r k ) = 



1 + e 



(1-e) 



(13) 



and averaging over different disorder realizations accord- 
ing to [l6| we have 



IV. TRANSMITTANCE AND QUENCHED 
LEVY WALKS 



Let us now turn to the quenched case. Within the 
scaling framework described by equations (|4I2I) . we first 
derive an analytic expression for the coefficients A e . We 
can exploit the fluctuation-response relation connecting 
Pa,e(r,t) to C Qj£ (L) i.e. the stationary conductivity of a 
system of size L. In particular according to [l9|, [22[ we 
have: 

C a , e {L)- 1 = lim / e iut (P a , s (L,t)-P a , e (0,t))dt (9) 

Plugging the scaling form © of P ay£ (r,t) into equation 
(|9]) and imposing that the scaling function depends on 
the transmittance e only through the constant A £ , as in 
(jlj, we obtain: 



C a , e (L) 



A\ +a L- a if0<a<l 
AlL~ x ifKa 



(10) 



where the proportionality constant is independent of e. 
Equation (fTU|) extends the Einstein relation for anoma- 
lous conductivity in 0, [22j, taking into account the 
transmittance e. The conductivity can be evaluated di- 
rectly by studying the stationary current in a system of 
size L and fixed boundary conditions. The (stationary) 



C a<e (L) 



D L+£ 
w n«(i- £ ) 

n 1+E 



if < a < 1 
if 1 < a 



(14) 



where the constant D a is independent of e. The results 
in equations (fT4| is consistent with (|4j and (fT0|) . with an 
analytical estimate for the multiplicative coefficients of 
the scaling lenghts: 



A F = 




if < a < 1 
if 1 < a 



(15) 



In Fig. [3] we compare our analytical prediction (1151) with 
numerical simulations at different value of a, with ex- 
cellent agreement. Clearly A e diverges for e — > 1, i.e. 
perfect transmission and A e vanishes for e — > — 1, i.e. 
total reflection. We remark that A e behaves differently 
than in the annealed model both in the superdiffusive 
and in the normal case. 



V. NUMERICAL EVIDENCES ON QUENCHED 
STRUCTURES 

Let us now present a numerical analysis of the 
quenched cases, evidencing that our assumption holds, 
i. e. that the whole effect of a variation in e can be sum- 
marized in a variation of the coefficient A £ , and therefore 
that the scaling functions are super-universal, i.e. f a {x) 
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FIG. 3: Value of the coefficient A E as a function of e. Symbols 
represent numerical simulations, while continuos lines are the 
analytical calculations of formula (fl5|) . 



in Eq. does not depend on e. For standard diffu- 
sion this properties obviously holds: indeed the scaling 
functions are Gaussian and they can be characterized ba- 
sically by their variance, i.e. the scaling length of the 
process. However in superdiffusive processes, for a < 1, 
f a (x) is a non trivial function decaying at large distances 
as x~ 1 ~ a , so an analogous property is not trivial. 

Numerical data, in Fig. [4] evidence for a = 0.7 in 
the random case that P Q , e (r, t) for different values of the 
transmittance can be scaled into a single function, inde- 
pendent of e. Moreover, the dashed line shows that the 
scaling function is different from the standard Levy func- 
tion describing the sum of independent Levy distributed 
random variables, even if they are characterized by the 
same long tail [f|. 

Let us now consider the random case with a > 1. The 
scaling length grows as in a diffusive process, and there- 
fore we expect the scaling function f a {x) to be a Gaus- 
sian, independently of e and a. Fig. [SJ obtained for 
a = 1.5, evidences this property for small r/£ e (t). How- 
ever, for large r/£ E (t), P a ^ E (r,t) contains now the sub- 
leading contribution h a _ £ {r,t). This term vanishes for 
large times but it can influence the momenta of the dis- 
tribution. Fig. [5]shows, at large r/£ E (t), the presence of 
this subleading term. In [18J, h a ,o(r,t) has been evalu- 
ated using a "single long jump" approximation. Here, we 
show that the same argument applies within the ansatz 
(fj| with the coefficients determined by (fTSjl . and leads 
to a correct estimate of the function h a ^ £ {r,t). In par- 
ticular, within the same approximation, for r/£ £ (t) 3> 1 
the probability of reaching a point at distance r is deter- 
mined by the probability of performing a single ballistic 
stretch of length r, times the number of scatterers visited 
by the walker in a time t; this number can be estimated 
as £ s (t)/A, where A is the average distance between the 
scatterers, that for a > 1 is finite and independent of e. 
If the whole effect of a variation in the transmittance e 



FIG. 4: Rescaling of P a ,e(r, t) for a — 0.7 on a random struc- 
ture. The value of £ e (t) are evaluated according to Q and the 
values of the coefficients A £ is given by (|15[). Red-dashed line 
represents the shape of the scaling function characterizing an 
uncorrelated Levy flight with the same a — 0.7. Notice that 
the tail of the distribution is the same since in both processes 
it is determined by the value of a. 

can be encoded in a change of A £ as in Eq. ([J}, then we 
expect that the behavior of h a ^{r, t) can be estimated as 

(16) 

where the proportionality constant is independent of e. 
In Fig. [6] we check indeed that our ansatz is correct and 
that the tail of the distribution P Q , £ (r, t) can be described 
according to equation (|T5|) . with the e-dependent coeffi- 
cient. 

Finally, we focus on the case of deterministic one- 
dimensional fractal quasi-lattices i.e. the lower panel of 
Fig. [TJ where the step length distribution is described 
by the parameters n u and n r . As explained in [l9j . for 
e = the motion of the random walker is ruled by the 
parameter a = log(n u )/log(n r ) which plays the role of 
a in the random structure. In particular, the scaling 
length of the process also in the deterministic case grows 
as for a < 1 and t x l 2 for a > 1. Here, we con- 

sider averages performed over processes starting from a 
given point of the structure, e.g the origin evidenced in 
Fig. [TJ For local quantities, the scaling function does not 
present long tails since arbitrary long jumps are placed 
far away from the starting point. On the other hand, for 
a < 1 the fractality of the structure induces character- 
istic log-periodic oscillations in the scaling function. In 
particular for e = we have [lj| 

P a ,o(r,t) = CW/4>/m.9(log„„4M)) (17) 

where £o(t) is the scaling length of the process on the 
quasi-lattice, f' a is the scaling function and g{x) is a 
function of period one. According to Eq. ((H) , a variation 
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FIG. 5: Rescaling of P a ,e(r, t) for a = 1.5 on a random struc- 
ture. The value of £ £ (t) are evaluated according to Q and the 
values of the coefficients A E is given by (|15[) . Dashed line rep- 
resents the gaussian behavior of the scaling function for not 
too large r/£ e (t). At larger r/£ e (t) there is a slowly decaying 
function h a ,e[r,t) superimposed to the gaussian behavior. At 
long times (black symbols) h a , e (r,t) tends to vanishes, how- 
ever it still influences the momenta of the distribution. 
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FIG. 7: Rescaling of P a ,e(r, t) for a — log(2) / log(3) on a frac- 
tal quasi-lattice. The value of £ £ (t) are evaluated according 
to (|4I15[) . Times are chosen according equation (fT9|) to detect 
scaling. 



log„ u £ £ {t) = k + log„ u £g(t) with k integer, i.e. 



(l-e)(l + e) 



(19) 



In Fig. [7] we evidence that for times chosen according 
to equation (TH)]) scaling holds for different e. The com- 
plex devil staircase shape of the scaling function is typical 
of these fractal structures. Notice that the choice of times 
flT9l) . as a function of e, is crucial to recover the scaling, a 
further test of the validity of Eq. (j4j . In the case a > 1 
the scaling functions are Gaussians, log-periodic oscilla- 
tion are absent and scaling can be recovered through Q 
and fll~5l) , without the tuning of times (jT9j , as in the ran- 
dom case. 



VI. TIME RESOLVED TRANSMISSION 



FIG. 6: Plot of -Pc*,£ {r, t) I £ £ (£) as a function of v for a = 1.5 on 
a random structure. The value of £ e (t) are evaluated accord- 
ing to Equation ([4]). The plot evidences that the tail h a ,e(r, t) 
of the distribution corresponds to Equation (|16|) where the ef- 
fect of different transmission coefficient can be summarized in 
the constant A E in (|15|l . 

of the transmittance only induces a rescaling of the cor- 
relation length, so that the scaling function for a generic 
e is expected to be : 

P«, e (r,i) = CH*)/a(r/4(«),Klog n „4(*))) (18) 

where £ e (t) and the corresponding coefficients are again 
given by equations (|4I15|) . 

Therefore for processes with different transmittances e 
and £, scaling holds if times t and t are chosen so that 



Due to its generality, the scaling approach can be use- 
ful not only in the analysis of P QjE (r, t) but also for other 
interesting physical quantities, such as the exit time prob- 
ability and the time resolved intensity fToj j , which are ex- 
perimentally more relevant. Let us consider the effect of 
a transmittance coefficient e ^ in the case of a ran- 
dom structure. We consider a walker starting at time 
t = from the border of a sample of size L and we 
define the time resolved transmitted intensity I a _ £ (t,L) 
as the probability for the walker to reach the bound- 
ary at distance L before returning to the starting point. 
The scaling hypothesis for the intensity I atE (t,L) reads 
I a , e (t, L) = Bg(tAl +a /L 1+a ) where g(-) is a scaling func- 
tion and the proportionality constant B depends on L, a 
and e. We calculate the coefficient B as follows. 

The conductivity C a<e (L) is by definition the total 
number of walkers escaping from a system of size L before 
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FIG. 8: The rescaled time-resolved transmitted intensity as 
a function of the rescaled time according equation (|2ip . on a 
random structure with a = 0.4. 

returning to the starting point, independently of time i.e. 

/>oo 

C a<s (L)= / I a>e (t,L)dt. (20) 
Jo 

Imposing the integral (|20|) to be proportional to A\ +a / L a 
according to equation (|10[) . we obtain 

J 0ie (t, L) = i- 2 "- 1 ^ 2 "^^^/^ 1 ^). (21) 

The scaling form (|2Tj) has been verified in Figure [8] for 
a = 0.4 and different values of e, with the values of the 
coefficients determined by (TPoj) . 

VII. CONCLUSIONS 



length (the parameter e in our model). This parameter 
enters in the multiplicative prefactor, ruling the depen- 
dence of the scaling length as a function of time. Estima- 
tion of the latter is of course relevant for finite samples 
and times. 

We obtained an analytic expression for the multiplica- 
tive coefficients of the scaling lengths, as a function of 
e in the annealed Levy walk case, evidencing its in- 
dependence on the transmittance in the superdiffusive 
regimes. Within the scaling framework, we also deter- 
mined a closed form for the coefficient in the quenched 
random and fractal cases, which are the most relevant for 
experiments. 

Another remarkable result of our analysis is that the 
scaling functions feature a superscaling property, namely 
they are independent of the transmittance and only de- 
pend on the structure through the exponent a. Interest- 
ingly, we employed the scaling properties to infer the de- 
pendence of the time-resolved transmission in finite sam- 
ples. All our analytic results have been compared with 
numerical simulations, with excellent agreement. 

As the scaling picture discussed here has been shown to 
hold also in higher dimensional cases [l(| , we expect that 
this superuniversality, holding in one dimension, can be 
detected also in higher dimensions. Moreover, we expect 
that superuniversality could hold not only for a varia- 
tion of the transmittance e but for a wider class of local 
transformation of the dynamics such as the introduction 
of waiting times or of second neighbors jumps. 

Besides their theoretical interest, our result are of im- 
portance to interpret correctly the experimental and nu- 
merical results. For instance, in the optical experiments 
of Ref. [3] it is possible to control the mean free path of 
light in diffusive media and investigate the approach to 
the scaling limits for the same distribution of glass sphere 
diameters. 



We have shown that Levy walk in a quenched, long- 
range correlated structures satisfies a generalized scaling 
relation for arbitrary values of the transmittance, both in 
the random and in the fractal expressed by Eq. 

© and (fTE)) respectively. The main difference between 
the two models is in the presence of the subleading term 
for the random case and of log-periodic oscillations in 
fractal quasi-lattices. As expected, all the leading scal- 
ing behavior are unaffected by a change of the scattering 
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